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ABSTRACT 

We report the discovery of rapid variations of a high-velocity C IV broad absorption line trough in 
the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the 
Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy 
were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant 
(> 4cr) variability in the equivalent width of the broad (~ 4000 km s -1 wide) C IV trough on rest-frame 
timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet 
reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two 
over the duration of the campaign. We evaluate several potential causes of the variability, concluding 
that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the 
outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale 
of variability places a lower limit on the density of the absorbing gas of n e ^ 3.9 x 10 5 cm -3 . 

The broad absorption line variability characteristics of this quasar are consistent with those observed 
in previous studies of quasars, indicating that such short-term variability may in fact be common 
and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy 
feedback scenarios. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — quasars: 
absorption lines 
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1. INTRODUCTION 

Broad absorption lines (BALs) in quasar spectra have 
observed velocity widths greater than 2000 km s _1 and 
are thought to be produced by high-velocity outflowing 
winds launche d near the q uasar’s su permassive black hole 
(SMBH; e.g., Weymann et al. 1991). Many theoretical 
investigations into these winds have supported a model 
in which the outflows are driven by radiation pressur e 


on lines (e.g., Murray et al. 1995 
although magnetoh ydrodynamical 


Proga et al. 2000) 

_'orces may contribute 

(e.g.,|Konigl & Kartje|1994|. These winds are potentially 
a means by which gas is evacuated from the host galaxy 
and thus have the potential to a ffect the growth of both 
the galaxy and its SMBH (e.g.,|Di Matteo et a l.| [2005 
Moll et al. 2007; Springel et al. 2005; King 2010). BALs 


arc observed in 10- 15%~o f optically selected quasars (e.g., 
Gibson et al.|2009[ Allen et al.|2011 ); as such, BALs play 
an important part in our understanding of both quasars 
and galaxy/quasar evolution. The detailed structures 
and locations of these outflows are not well determined, 
and we must understand their environments to deter¬ 
mine if BALs are viable agents of feedback in their host 
galaxies. 

Variability of BAL troughs provides a powerful tool 
with which to constrain the properties of these out¬ 
flows. BALs are observed to be variable in time, both in 
strength (frequently characterized by equivalent width) 
and in shape. There are many possible causes of the vari¬ 
ability, including movement of the absorbing gas (“cloud 
crossing”) or changes in the ionizing radiation received 
by the gas. Several studies show support for the vari- 
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ability ar ising; due to bulk movement of the outflowing 


gas (e.g., Lundgren et al. 2007 

Gibson et al. 

2008 Hall 

et al. 

2011 

Vivek et al. 2012 

Capellupo et al. 

2013), 


changes or changes in “shielding gas” as the cause (e.g., 
the coordinated variabilit y seen in multiple troughs of 


the same transition by Filiz Ak et al. 2012 2013). It 


is possible that the mechanisms causing variability dif¬ 
fer from object to object, yet each mechanism allows one 
to constrain different properties of the outflow environ¬ 
ment. For example, if the variability arises from cloud 
crossing, one can estimate the radius of the absorbing 
gas from the ionizing continuum source. If the variabil¬ 
ity is due to ionization response, one can estimate the 
density of the outflow. Either of these parameters can 
help determine whether or not BALs are viable sources 
of large-scale feedback to their host galaxies. 

The variability timescales probed by previous BAL 
studies vary widely — studies have repor ted BA L-trough 
varia bility on timescales of years (e.g., 

2013|) all the way down to timescales 


Filiz Ak et 
of 8-10 


al. 

days 


(|Capellupo et ah 2013), with many studies examining 


times cales in between (e.g.,|Barlow|1993; Lundgrcn et al. 
2007|. BAL variability has been seen on all of the 


seen on 

timescales probed; however, to date, no studies have re¬ 
ported variability on rest-frame timescales shorter than 


be less variable on shorter timescales (e.g., 

Gibson et al. 

2010 Filiz Ak et al. 

2013). Detection of variability on 


ability mechanisms and outflow environments. 

Recent spectroscopy from the Sloan Digital Sky Survey 


Reverberation Mapping Project (SDSS-RM; Shen et al. 


20151, whose primary goal is reverberation mapping of 


the broad emission line regions of quasars, also includes 
observations of a number of quasars hosting BALs. The 
cadence of these observations allows investigation of BAL 
variability in a number of targets down to short rest- 
frame timescales. In particular, preliminary measure¬ 
ments for one target indicated Civ BAL variability on 
short (~1 day) rest-franre timescales. This target, the 
quasar SPS S J141007.74+541203.3 (also known as SBS 
1408+544A; Chavushyan et al.|1995|), was o b served spec¬ 
troscopically liTTMTTjrSfepamanetaL (1998), who 
noted strong absorption in the C IV emission-line region. 
There was also one previous spectrum taken as a part 
of the S DSS-III Baryon Oscillation Spectroscopic S urvey 
(BOSS; |Eisenstein et al.||2011| |Dawson et al.||2013|) from 
2013 M ay that was releas ed in the SPSS Data Release 12 
(DR12; Alam et al.||2015l). Based on the Mgn and Chi] 
emission lines from the PR12 spectrum, the redshift of 
this quasar is z = 2.337 ± 0.00 3. It has an apparent 
i-band magnitude m* = 18.1 (Al am et ah1[2015) and ab¬ 
solute magnitude Mi = —27.69 fPans et ai.|[2UI5). 

We here report on the detailed investigation of this 
target and the implications of these measurements for 
the environment producing the BAL trough. In Sec¬ 
tion^ we present the data used in our investigation and 
its preparation for analysis. We discuss our measure¬ 
ments of the BAL parameters and the search for vari¬ 
ability in Section [3] and follow with a review of possi¬ 
ble causes of variability and the physical constraints we 
can derive from our data in Section ^ We conclude in 
Section [5] by summarizing our results and their impli¬ 


cations. Throughout this work, we assume a cosmol¬ 
ogy with H 0 = 70 km s -1 Mpc -1 , 11 m = 0.3, and 
P A = 0.7. 

2. DATA AND DATA PREPARATION 
2.1. BOSS Spectra 

The spectra utilized in our study were acquired as 
a part of the SDSS-RM project, which is a dedicated 
multi-object reverberation mapping campaign performed 
as a part of SDSS-III BOSS. For a tec hnical overview 
of the SDSS-RM project, 


see 


Shen et al. (2015). In 


brief, 849 quasars were spectroscopically monitored with 
the BOSS spectrog r aph on the SPSS 2 5-meter telescope 


(Smee et al. 2013 
uary 1 


Gunn et al. 2006) from 2014 Jan- 
through 2l)l'4 July 3, resulting in 32 epochs of 
spectroscopic observations over this period. The median 
spacing between observatio ns is about 4 days — see Ta¬ 
ble 2 of Shen et al.| (2015) for a complete log of these 


observations. The 10 V of the BOSS spectrograph is 2.5 
degrees in diameter, and its wavelength coverage is 3650- 
10,400 A with a spectral resolution of R ~ 2000. Because 
of the need for high-precision relative spectrophotometric 
calibrations for the SDSS-RM project, the data were pro¬ 
cessed using an improved flux-calibration method that 
uses extra simul taneous sta ndard star observations (see 
§3.3 and 3.4 of Shen et al. 20151. The first SDSS-RM 
BOSS spectrum was obtained on MJD 56660.209 — we 
define this as the “beginning” of our campaign, and all 
light curves presented as a part of our analysis (unless 
explicitly stated otherwise) show the rest-frame time be¬ 
tween each epoch and MJD 56660.209 (2 January 2014). 

The processed BOSS spectra from this campaign were 
first corrected for Galactic ext inction using a Ry = 3.1 
Milky Way exti nction model (Carde lli et al. 1989) and 
Ay values from Schlegel et al. (1998). Before beginning 
analysis, we converted the observed wavelengths of the 
spectra to the rest frame. All 32 epochs of data were 
visually inspected for quality; the median signal-to-noise 
ratio (S/N) of the spectra was ~30 per pixel, with most 
spectra falling between 20-40. We excluded one spec¬ 
trum from our analysis due to low S/N (< 10), leaving 
31 SDSS-RM spectral epochs for our analysis and the 
earlier observation from the DR12 data release. Figurejl] 
shows the mean spectrum of our quasar duri ng the S DSI+ 
RM campaign (RMID 613 in Table 1 of She n et al.|2015 ). 
Visual inspection of the spectrum reveals” BAL quasar 
with low reddening and a prominent C iv BAL trough 
that is at sufficiently high velocity to be detached from 
the CiV emission line. There are no BAL features in 
lower-ionization transitions like Mg II or Al ill, making 
this target a high-ionization BAL quasar. There are also 
two narrower C iv absorption features at lower outflow 
velocities superimposed on the C IV emission line, as well 
as N V absorption redward of the Lyman-a emission line, 
but this region is heavily contaminated by intervening 
Lyman-a absorption and thus its analysis is outside the 
scope of this work. There is a hint of broad Si IV ab¬ 
sorption at a similar velocity as the C IV BAL, but it is 
heavily contaminated by Cn and Oi/Sin features and 
does not ever reach a depth below 90% of the continuum 
level. 

We also inspected the spectrum taken on 1991 October 
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Fig. 1.— Left panel: The mean spectrum (black solid line) from all 31 epochs used in this analysis with its continuum fit (blue dashed 
line). The line-free regions chosen to include in the continuum fit are shaded in gray and also highlighted in red. The approximate 
wavelength coverage of the SDSS g- and z-band photometric filters are shown as solid b lu e horiz ontal lines. Right panel: The SDSS-RM 
spectrum (black) from MJD 56683 and its accompanying PCA fit (blue) from|Paris et al.|(|2015|). We use the PCA fit to remove the Civ 
emission-line signatures, and investigate the lower-velocity C IV absorption features. 



Fig. 2.— The 1991 spectrum from |Chavushyan et al.| (1995), 
smoothed by 5 pixels and normalized by the mean tlux density 
(black solid line) and our mean continuum-normalized spectrum 
from the SDSS-RM campaign (red solid line). The blue dashed 
lines show the rest-frame wavelengths of N v, Si IV, and C IV for 
reference. The red dotted lines show the wavelength region spanned 
by the C IV BAL trough during the SDSS-RM campaign. 


4 by Chavushyan et al. (1995) for comparison (Figure[2|). 
The (JIV UAL trough does not appear to be present m 
the 1991 spectrum, though the S/N is too low to de¬ 
termine definitively whether or not the BAL trough is 
present; however, the depth measured within the trough 
region in the 1991 spectrum is consistent with there be¬ 
ing no trough flux below 90% of the continuum level to 
~3cr significance. The lowest-velocity, narrower C IV ab¬ 
sorption feature close to the ClV fine center is clearly 
present in the 1991 spectrum, though we cannot tell if 
the nearby higher-velocity narrow absorption feature su¬ 
perimposed on the C iv emission fine is present due to 
low S/N. 


2.2. Continuum Fits and Normalization 

To disentangle the BAL variability from the overall 
quasar variability, we fit a continuum to each individual 


spectrum. Following §3.1 of Filiz Ak et al. (20121, we 
first attempted to fit a reddened power law model to the 
continuu m emis sion using the SMC-like reddening model 
from Pei (1992) - however, the results were consistent 

with no extinction, so a simple, non-reddened power law 
yielded suitable fits for this particular target. We used a 
nonlinear least-squares algorithm and fit to four regions 
within the spectra that were relatively free of emission 
and absorption features. We chose these relatively line- 
free regions, shown in Figure[l] by visually inspecting the 
mean spectrum an d comparing it wit h the SDSS compos¬ 
ite spectrum of|Vanden Berk et al.|([2001). We used the 
same line-free regions for each epoch. Because some of 
the line-free regions contained more pixels than others, 
we weighted each pixel such that each individual region 
contributed equally in the continuum fit. This ensures 
that none of the line-free regions holds more weight in 
the fit just because it spans a larger wavelength range 
(and thus contains more pixels). 

To characterize the uncertainties in the continuum 
fits, we account for two different sources of error. The 
first component accounts for the uncertainty caused by 
noise in the spectra. To evaluate these uncertainties, we 
used “flux randomizati on” Monte Carlo iterations as de¬ 
scribed by Peterson et al. (19981), where the flux in each 


pixel of the spectrum is altered by a random Gaussian 
deviate based on the uncertainty assigned to that partic¬ 
ular pixel. The continuum is then fit to the new spec¬ 
trum, and the process is repeated 100 times. We adopt 
the standard deviation of these 100 iterations as the first 
component of the uncertainties in the continuum fit. 

The second component to the continuum uncertain¬ 
ties considers the small changes in the fit resulting from 
the decisions we made when fitting the continuum. For 
example, we choose which particular line-free regions of 
the spectrum are used in the fit, and whether or not 
each pixel is given equal weight. To evaluate the un¬ 
certainties introduced by these decisions, we created two 
more continuum fits for each spectrum. First, we fit the 
continuum in each epoch using two additional line-free 
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Fig. 3. — Mean and RMS residual spectra for the continuum-normalized (left panels) and of the PCA-normalized (right panels) cases. 
The PCA fit generally results in similar mean and residual spectra compared to our continuum fit, but also allows one to isolate and 
characterize the low-velocity C IV absorption near the broad emission. Each panel shows a different wavelength region along the spectrum. 
The top subpanel for each region displays the mean continuum-normalized spectrum over the 31 epochs, and the bottom panel shows the 
RMS residual spectrum. Dotted horizontal lines across the mean spectrum show a normalized flux density of 1.0 and 0.9 to help the eye in 
identifying absorption features. Vertical dashed lines indicate the wavelength regions spanned by the Civ BAL (middle panels) and also 
in which absorption associated with the Civ BAL would be expected to appear for Nv and Si IV (top panels) and Mg II (bottom panels). 


regions, but still using the same weighting scheme as de¬ 
scribed above, where each line-free region is given equal 
weight. Second, we fit a continuum using our original 
line-free regions while weighting each pixel equally. We 
then adopted the average deviation of these two contin¬ 
uum fits from the original continuum fit as the second 
component to the continuum uncertainty. 

The two uncertainty components are added in quadra¬ 
ture and adopted as our formal uncertainties in each con¬ 
tinuum fit. Although a continuum was fit to each epoch 
individually, Figure [l] shows the fit to the mean spectrum 
as an example of how well the fits performed. Each indi¬ 
vidual epoch was divided by its continuum fit to obtain 
a normalized spectrum, and both the spectral and con¬ 
tinuum fit uncertainties were propagated to determine 
the final uncertainty on the normalized spectrum. Here¬ 
after, we will refer to these normalized spectra as the 
“continuum-normalized” spectra. In all cases, the uncer¬ 


tainties in the continuum-normalized spectra are domi¬ 
nated by the spectral uncertainties rather than the for¬ 
mal uncertainties in the continuum model. We created 
a mean and root mean square (RMS) residual spectrum 
from the continuum-normalized spectra; the RMS resid¬ 
ual spectrum is obtained by taking the standard devia¬ 
tion about the mean spectrum, the result of which shows 
the variable components (Figure^]). The base level of the 
RMS flux (on the order of ~ 0.2) in all panels is above 
zero due to noise in the spectra, which adds a constant to 
the entire RMS spectrum. The RMS spectrum contained 
a few sharp features arising from poor/inconsistent sky 
subtraction; these small segments of the spectra were 
masked to avoid confusion, as they did not lie near any 
of the emission or absorption features of interest in this 
study (all fell between 1600-1900 A in the rest frame). 
The region of the spectrum redward of ~2200 A suf¬ 
fered from telluric contamination; these features were too 
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widespread to mask, so they are shown in full in Figure[3} 


2.3. Principal Component Analysis Fits 

As a part o f the SDSS DR12 quasar catalog produc¬ 
tion (Paris et al. 2015), a principal component analysis 
(PCA) was done on each epoch of BOSS spectroscopy. 
This procedure was also used for the DR9 and DR10 
quasar catalogs (Paris et al.|2012 2014); see these works 
for a description ot this technique! PCA fits both the 
emission-line and continuum features in each spectrum 
using a reference sample of quasars for which the rele¬ 
vant features are observed. An example of a PCA fit to 
one of our spectra is shown in Figure [T] In this target, 
two lower-velocity and narrower C IV absorption features 
are seen superimposed on the C IV emission line, and the 
C IV emission-line flux in this low-velocity region is well- 
fit by the PCA approach. Thus, we consider the PCA 
fits in our analysis to help eliminate contributions from 
the C iv emission line, particularly relevant for measure¬ 
ments of the lower-velocity absorption features. As with 
our continuum fits, we divide the original spectra by their 
PCA fits to obtain the “PCA-normalized” spectra at 
each epoch for use when investigating the lower-velocity 
C IV absorption features. Figure [51 shows the mean and 
RMS residual spectrum of the PCA-normalized spectra. 
All of the emission features are well-removed from the 
mean and RMS residual spectra by the normalization 
except for the Lyman-a region, which is blended with 
N V and has significant Lyman-a absorption on top of 
the emission line. The poor fit is due to the sample size 
used in the PCA in this region below rest-frame 1280 
A, the sample used to perform PCA is small and makes 
use of hand-fitted continua, so we are ab le to repr oduce 
fewer features in this region (e.g., |Paris et al. 2012]) . As 
with the continuum-normalized spectra, we mask a few 
regions with significant sky features in the RMS residual 
spectrum to avoid confusion. 

2.4. Photometric Observations 

As a part of the SDSS-RM program, photometric ob¬ 
servations of our target were also obtained with the 
2.3m Steward Observatory Bok telescope and the 3.6m 
Canada-France-Hawaii Telesco pe (CFHT) . The observa¬ 
tions are described in §3.5 of Shen et al. (20151. The 
Bok/90Prime instrument has a 1 deg 2 field of view and 
a plate scale of 0745 pixel -1 . About 60 nights were allo¬ 
cated to the SDSS-RM program from January-July 2014, 
yielding 31 epochs in the g band and 27 in the i band by 
the end of the program. The CFHT observations were 
taken with MegaCam, which also has a 1 deg 2 field of 
view, with a pixel size of 07187. Over the course of the 
campaign, 26 epochs were observed in the g band and 
20 observe d in the i band with the CFHT (see |Fukugita 
et al. 1996 for a description of the ugriz filtersjl 

In this target, the a filter covers both Lyman-a/Nv 
and ClV (see Figure [l]), both of which are observed to 
vary during the campaign, particularly the Lyman-a/N V 
region (Figured). Thus, we do not consider the g-band 
photometric lignt curve to be a good proxy for the contin¬ 
uum variability. However, the i-band photometric light 
curve covers a relatively line-free region of the spectrum 
for this object (Figure TT1) so we include this light curve 
in our investigations (Figure [4]) . Aperture photometry 
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Fig. 4. — Photometric i-band continuum light curves from the 
CHFT (black filled circles) and Bok (red squares) telescopes. The 
BAL quasar varies by only ~10% over this time period. 


was performed o n all images using SExtractor (jBertin 
fe Arnouts||1996) and a preliminary flux calibration was 


performed. A zero-point shift was applied to place the 
Bok photometry onto the CFHT magnitude scale, which 
in turn is tied to the photometry of SDSS objects in our 


field a s a step in the MegaPipe reduction pipeline (Gwyn 


20081. Several epochs were identified as outliers for ail 


targets observed and thus categorized as bad data and re¬ 
moved from the light curves; we had good measurements 
for 43 total epochs. The quasar did not exhibit very 
much photometric variability during the campaign; the 
flux changed by a factor of only about 10-15% during 
this period. Measurements of the flux densities within 
the line-free regions of the spectra similarly yield light 
curves consistent with little variability. 


3. VARIABILITY MEASUREMENTS 
3.1. Searching for Variability 

Inspection of the C IV emission-line region in the mean 
and RMS residual of both the continuum-normalized and 
PCA-normalized spectra (Figures [3] and [5|) reveals sig¬ 
nificant variability in the high-velocity ClV absorption 
trough centered near —16700 km s -1 (hereafter referred 
to as Trough A; see Figure [5b. Figures [3] and [5] also 
show some variability in the (J iv emission line, which 
makes it more difficult to determine whether there is any 
variability in the lower-velocity C iv absorption troughs 
using the continuum-normalized spectra. However, the 
RMS spectrum of the PCA-normalized spectra largely 
removes the variability signatures from the emission-line 
features (Figure [3], allowing us to investigate variabil¬ 
ity in the lower-velocity absorption features. There is 
detectable low-amplitude variability in the associated 
medium-velocity Civ trough (hereafter Trough B; cen¬ 
tered near —4800 km s _1 ), but none is detected in the 
lowest-velocity associated narrow C IV absorption trough 
(hereafter Trough C; centered near —1000 km s _1 ) su¬ 
perimposed on the ClV emission line. The RMS resid¬ 
ual signal in the Trough A region is nearly identical to 
that in the continuum-normalized RMS spectrum, which 
is expected given that this trough is detached from the 
C IV emission line and thus we expect little contamina¬ 
tion from the emission-line variability. 

Examining the RMS residual spectra derived from 
both the continuum-normalized spectra and the PCA- 
normalized spectra at various other wavelength regimes 
allows a search for associated (or non-associated) absorp¬ 
tion and emission-line variability across the spectrum. 
There are signs of variable absorption in the N V, Al ill, 
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Wavelength (A) 


Fig. 5.— Mean and RMS residual spectra for the continuum- 
normalized case, expanded to emphasize the Civ region of the 
spectrum. Dotted horizontal lines across the mean spectrum show 
a normalized flux density of 1.0 and 0.9 to help the eye in identi¬ 
fying absorption features. 

and Si IV regions (as evidenced by a slight increases in 
the RMS residual flux in these regions), but these sig¬ 
natures are not significant enough to disentangle them 
properly from the noise, and as discussed in Section [2.2[ 
none hosts a BAL trough of >10% fractional depth. We 
hereafter focus our attention on the ClV region of the 
spectrum shown in Figure [5] 

3.2. Measurements 
3.2.1. BAL-Trough Properties 

To determine the minimum and maximum veloci¬ 
ties of Trough A, we smoothed the mean continuum- 
normalized spectrum using a boxcar smoothing algo¬ 
rithm over five pixels. As is standard in BAL stud¬ 
ies, we then determined the velocities between which 
the flux in Trough A remained below 90% of the 
continuum level (i.e., where the normalized flux den¬ 
sity is < 0.9) and obtained u m i n = —14470 km s -1 and 
Umax = —18810 km s-\ which yields a trough velocity 
width of 4340 km s- . The minimum and maximum ve¬ 
locities of the trough in each epoch were measured sepa¬ 
rately, and were consistent with those measured for the 
mean spectrum to within the uncertainties (the range 
of measured values was within ~200 km s -1 of the 
mean values). We applied the velocity limits found in 
the mean spectrum to measure the rest-frame equivalent 
width (EW), mean fractional depth, and absorbed-flux- 
weighted centroid velocity (u cen t) of Trough A at each 
epoch using the non-smoothed continuum-normalized 
spectra (Table [I]). Uncertainties in each parameter were 
calculated using flux-rando miza tion Monte Carlo itera¬ 
tions as described in Section [272] We adopt the standard 
deviation of these iterations as our measurement uncer¬ 
tainties and tabulate them in Table |T] 

Figure [g] shows the EW, mean fractional depth, and 
lucent measurements of Trough A as a function of rest- 
frame time from the beginning of the RM campaign. 
Both the trough EW and depth increase from At = 0 
to ^15 days by about 23%, followed by a steady decline 
over the rest of the campaign by a factor of ~2. The 
similarity of the EW and mean fractional depth light 
curves indicates that the trough profile does not change 


TABLE 1 
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MJD 

Af rest b 
(days) 

EW 

(A) 

Mean 

Depth 

^cent 

(km s -1 ) 

56426.00 

0.00 

4.37 

± 

0.08 

0.203 

± 

0.004 

—16799± 25 

56660.21 

70.19 

5.07 

± 

0.08 

0.235 

± 

0.004 

— 16768=t 20 

56664.51 

1.29 

4.80 

± 

0.11 

0.223 

± 

0.005 

—16768± 30 

56669.50 

1.49 

4.77 

± 

0.16 

0.221 

± 

0.007 

— 16698=t 41 

56683.48 

4.19 

5.06 

± 

0.08 

0.234 

± 

0.004 

—16695± 17 

56686.47 

0.90 

5.37 

± 

0.08 

0.249 

± 

0.004 

—16707± 20 

56696.78 

3.09 

5.78 

± 

0.08 

0.268 

± 

0.004 

—16707± 17 

56715.39 

5.58 

6.24 

± 

0.07 

0.289 

± 

0.003 

—16706± 17 

56717.33 

0.58 

5.81 

± 

0.08 

0.269 

± 

0.003 

—16711± 15 

56720.45 

0.93 

6.02 

± 

0.09 

0.279 

± 

0.004 

—16724± 17 

56722.39 

0.58 

5.95 

± 

0.06 

0.276 

± 

0.003 

—16737± 11 

56726.46 

1.22 

5.81 

± 

0.09 

0.269 

± 

0.004 

—16727± 16 

56739.41 

3.88 

5.64 

± 

0.08 

0.261 

± 

0.004 

— 16688=t 19 

56745.28 

1.76 

5.28 

± 

0.08 

0.245 

± 

0.004 

—16687± 19 

56747.42 

0.64 

5.09 

± 

0.08 

0.236 

± 

0.004 

-16743± 25 

56749.37 

0.59 

5.40 

± 

0.09 

0.250 

± 

0.004 

— 16721=t 19 

56751.34 

0.59 

5.49 

± 

0.08 

0.254 

± 

0.003 

—16735± 18 

56755.34 

1.20 

4.95 

± 

0.09 

0.229 

± 

0.004 

—16741± 21 

56768.23 

3.86 

5.22 

± 

0.07 

0.242 

± 

0.003 

— 16756=t 17 

56772.23 

1.20 

4.79 

± 

0.07 

0.222 

± 

0.003 

-16729± 19 

56780.23 

2.40 

4.41 

± 

0.09 

0.205 

± 

0.004 

-16746± 22 

56782.25 

0.60 

4.32 

± 

0.08 

0.200 

± 

0.004 

— 16768=t 22 

56783.25 

0.30 

4.51 

± 

0.08 

0.209 

± 

0.004 

—16755± 21 

56795.18 

3.57 

4.01 

± 

0.09 

0.186 

± 

0.004 

— 16776=t 25 

56799.21 

1.21 

4.07 

± 

0.08 

0.188 

± 

0.004 

—16783± 22 

56804.19 

1.49 

3.90 

± 

0.08 

0.181 

± 

0.004 

—16774± 25 

56808.26 

1.22 

4.08 

± 

0.08 

0.189 

± 

0.003 

— 16739=t 25 

56813.23 

1.49 

3.77 

± 

0.08 

0.175 

± 

0.004 

—16750± 27 

56825.19 

3.58 

3.68 

± 

0.07 

0.171 

± 

0.003 

—16790± 29 

56829.21 

1.21 

3.34 

± 

0.08 

0.155 

± 

0.004 

— 16866=t 33 

56833.21 

1.20 

3.41 

± 

0.08 

0.158 

± 

0.004 

—16782± 32 

56837.19 

1.19 

3.65 

± 

0.07 

0.169 

± 

0.003 

— 16796=t 24 

mean 0 


4.83 

± 

0.01 

0.224 

± 

0.001 

— 16739=t 4 


8 All measurements are made with the continuum-normalized spec¬ 
tra. 

° All At measurements are made with respect to the previous 
epoch in the rest-frame. 

c Measurements made from the mean spectrum created from the 
31 SDSS-RM spectra — excluding the early spectrum from DR12. 

significantly over the period of the campaign; we investi¬ 
gate this further below. The centroid velocity decreases 
slightly over the campaign, but with a small amplitude 
of ~ 100 km s -1 that corresponds to less than two pixels 
in the spectral direction. 

We also measure the EW in all three C IV absorp¬ 
tion troughs (A, B, and C) using the PCA-normalized 
spectra, as the latter two cannot be examined with the 
continuum-normalized spectra due to their positions on 
top of the C IV emission line. We examine the result¬ 
ing light curves (Figure [7]): The widths of Trough B 
and Trough C are ~1311 km s -1 and ~1380 km s -1 , 
respectively, classifying them as mini-BALs. For our pur¬ 
poses, we do not need to consider each component of the 
ClV doublet separately in these troughs. Nonetheless, 
we note that Trough B appears to be a line-locked sys¬ 
tem of two C IV absorbers, with the weaker one closer 
to the systemic redshift. The light curve for Trough A 
created from the PCA-normalized spectra is very similar 
to that obtained using the continuum-normalized spec¬ 
tra; there is again a clear increase, followed by a long¬ 
term decrease in EW. However, the PCA fits produce 
more scatter in the EW light curve of Trough A than 
the continuum fits, likely due to uncertainty in deter¬ 
mining the boundaries of the BAL trough for exclusion 
during the PCA fitting. Thus, we rely on our continuum- 
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Rapid C iv BAL Variability 



Rest-frame days since MJD 56660.209 

Fig. 6.— Equivalent width (top panel), mean fractional depth 
(middle panel) and centroid velocity (bottom panel) of the high- 
velocity Civ BAL trough (Trough A; centered near —16700 
km s -1 ) as a function of days since the RM campaign began, in 
the rest frame. The point at far left is the additional observation 
of our target taken by SDSS, about 70 rest-frame days before the 
SDSS-RM project began. The apparent decline in velocity cen¬ 
troid, on the order of ~100 km s — , is of much smaller amplitude 
than the changes seen in the EW and depth, and corresponds to 
less than two pixels in the spectral direction. 



0 10 20 30 40 50 

Rest-frame days since MJD 56660.209 


Fig. 7. — Equivalent widths of Trough A (centered at ~ —16700 
km s _1 ), Trough B (centered at ~ —4800 km s' 1 ), and Trough C 
(centered at ~ —1025 km s' 1 ) absorption components as a function 
of rest-frame days since the RM campaign began. These measure¬ 
ments were made using the PCA-normalized spectra rather than 
the continuum-normalized spectra, as Troughs B and C lie on top 
of the C IV emission line. 


normalized spectra for our quantitative assessment of 
Trough A. Our conclusions below about the BAL vari¬ 
ability still hold when considering the PCA-normalized 
spectra — however, we opt for the conservative route of 
using the continuum-normalized spectra for our quanti¬ 
tative analysis. However, when comparing the behavior 
of Trough A with Troughs B and C, we use the PC A fits 
for all three troughs to avoid introducing possible sys¬ 
tematic uncertainties between the continuum-normalized 
and PCA-normalized spectra. 

As expected based on the RMS of the PCA-normalized 
spectra, there is no significant variability of Trough C, 
but the EW light curve for Trough B shows a similar be¬ 
havior to that of the light curve of Trough A (though with 
smaller fractional variability). The EW light curves for 
Trough A and B from the PCA-normalized spectra have 
a Spearman rank correlation coefficient of 0.525 with a 
p-value of 0.002, indicating that the two are correlated. 
The cross-correlation function for the two light curves 
(Figure |8]), showing the Pearson coefficient as a func¬ 
tion of time delay between the two light curves, similarly 
reveals that they are well-correlated, with a peak r ~ 
0.7. We are unable to measure any time delay between 
the two troughs with sufficient precision to be useful; we 
measure a time delay of 17 ± 19 days, which is consistent 
with zero. 


3.2.2. Short-Term EW Variability 

To search for short-timescale variability between 
epochs in Trough A, we calculated the change in EW 
(AEW) between each pair of sequential observations us¬ 
ing the continuum-normalized spectra. Uncertainties in 



Rest-frame Time Delay (days) 


Fig. 8. — The cross correlation function (CCF; black solid line) 
between the EW light curve of Trough A and the EW light curve of 
Trough B, showing the Pearson’s correlation coefficient as a func¬ 
tion of time delay. The auto-correlation function (ACF), which 
shows the high-velocity trough EW light curve correlated with it¬ 
self, is show for Trough A (red dashed line) and Trough B (blue 
dotted line). To eliminate contributions from the Civ emission 
line, these calculations were done using light curves made from the 
PCA-normalized spectra. 
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AEW (<taew) were calculated by summing the uncer¬ 
tainties of the two co ntributing EW measur ements in 
quadrature. Following Capellupo et al. (20131, in order 
to be considered “real” variability, we require a signifi¬ 
cance of 4 a (that AEW is at least four times as large 
as ctaew)- We identified four instances of secure vari¬ 
ability between sequential epochs during the SDSS-RM 
campaign, with variability on timescales down to 1.20 
days (~ 29 hours) in the quasar rest frame. Figure [9] 
displays the four variable pairs of spectra. 

To assess if this short-term variability in the trough is 
consistent with coordinated variations across the entire 
trough, we took the ratio between each pair of spectra 
and determined whether the flux ratio in the Trough A 
region is consistent with a constant to within the mea¬ 
surement uncertainties. The flux ratios are shown in Fig¬ 
ure [lO] along with a constant fit to the region across the 
trough. In all cases, a constant model fits the flux ratio 
in the trough region quite well. However, to be sure that 
the two sub-troughs within Trough A are not varying dif¬ 
ferently, we measured individual EW light curves for the 
blue half and red half of Trough A (which can be decom¬ 
posed into two BALs of widths at least ~ 2470 km s -1 
and ~1870 km s^ 1 , respectively, when ignoring the over¬ 
lapping wings). We find that the two sub-troughs show 
the same behavior; i.e., when the blue half is increasing in 
EW, so is the red half, and by similar fractional changes 
in EW. Thus, the trough is varying in a very similar, or 
coordinated, fashion acro ss all v elocities. 

As discussed in Section [3. 2. 1[ examination of the EW 
light curve of Trough A (Figure [6]) shows a steady in¬ 
crease in EW from rest-frame days 0-15, followed by a 
steady decreasing trend for the remainder of the cam¬ 
paign. This behavior raises the concern that our ob¬ 
served short-term variability could be merely caused by 
statistical noise about this long-term trend, which would 
affect the interpretation of the variability. Three of the 
four pairs of variable epochs are found in the decreasing 
half of the light curve, so to determine whether or not 
these features are likely to be statistical noise, we fit a 
linear trend to the light curve between days 15-55. A 
simple linear fit did not perform particularly well (the 
fit has a reduced y 2 of 5.02), so we instead fit the linear 
trend allowing for an intrinsic scatter com ponen t. We use 
the Bayesian linear-regression model of Kelly (2007) to 
fit the line, and find an intrinsic scatter of 0.18 A, which 
is twice the typical measured uncertainty in EW. Ex¬ 
plaining this scatter without intrinsic variability would 
require that our EW uncertainties are underestimated 
by at least a factor of two: an error this size is unlikely 
because we account for both pixel uncertainties and noise 
intr oduc ed by the continuum-fitting procedure (see Sec¬ 
tion 2.21. We present the linear fit, intrinsic scatter, and 
uncertainties in the fit in Figure ED 

To explore further the likelihood that our observed 
short-term variability is due to genuine variability, we 
created a simulated EW light curve for Trough A, plac¬ 
ing each epoch directly on the fitted linear trend and 
assigning uncertainties equal to our measured uncertain¬ 
ties at each epoch. We then allowed the EW in each 
epoch to deviate by a random Gaussian deviate reflect¬ 
ing these uncertainties, and counted the number of epoch 
pairs that deviated from each other by greater than 4cr 


solely due to statistical fluctuations. This simulation was 
designed to mimic our actual search for pairs of epochs 
that varied in our data. We repeated this test 10 6 times 
to determine the probability that the observed varia¬ 
tions are due simply to statistical fluctuations about the 
long-term decreasing trend. The simulations indicate a 
probability of only 0.025% for statistical noise to pro¬ 
duce three pairs of sequential EWs that differ by at least 
4 <t. When we restrict the simulation to count only those 
pairs with A t < 2 days in the rest-franre (of which we 
found two pairs in our light curve), the probability drops 
to 0.0076%. Thus, it is highly unlikely that the detected 
variations are due to statistical fluctuations about the 
long-term trend, and we conclude that they represent 
real short-term intrinsic variability in this system. 

3.2.3. Trough-Profile Variability 

Though the variations in trough-profile shape are min¬ 
imal on the shortest timescales, we also search for indi¬ 
cations of variability in trough-profile shape associated 
with the underlying longer-term trends seen in EW in 
Trough A. To do this, we divided spectra into groups con¬ 
taining spectra that were close together in time and had 
similar EW measurements. We averaged the continuum- 
normalized spectra in each group to improve the S/N, al¬ 
lowing us to compare the average trough profiles at differ¬ 
ent times throughout the campaign. Grouping the spec¬ 
tra resulted in six groups, each containing 3-5 sequential 
spectra with similar EW values. Group 1 contains the 
spectra taken from MJD 56660-56669, Group 2 between 
MJD 56717 -56722, Group 3 between MJD 56745-56752, 
Group 4 between MJD 56780-56783, Group 5 between 
MJD 56795-56808, and Group 6 between MJD 56813- 
56837. We averaged all of the normalized spectra in 
each group and plot each mean spectrum in Figure |12| 
The average trough profiles of each group appear simi¬ 
lar throughout the campaign; it is only the strength of 
the trough that appears to change. To test this conclu¬ 
sion, we measured the flux ratios between each pair of 
group-averaged spectra at each pixel of the trough. We 
find that in all cases, a constant model fit to the flux ra¬ 
tio between the different group-averaged spectra is again 
a statistically acceptable fit, indicating that the absorp¬ 
tion trough is changing in a coordinated manner across 
its entire velocity span on longer timescales as well. 


3.2.4. Emission-Line Variability 

The SDSS-RM campaign also offers the opportunity to 
search for links between emission-line and BAL-trough 
variability and thus possible connections between the 
high-ionization broad-line region and BAL outflows; such 
connections have been suggested by many authors (e.g., 
Chiang fe Murray 19961 Richards et al. 2011[ Hig-| 
ginbottom et al.| 2013|). As such, we also measured 
the integrated flux within the C iv emission li ne d ur- 
ing the SDSS-RM campaign, shown in Figure [l3j to 
look for possible coordinated variability. We calculated 
the flux in only the red half of the C IV emission line 
(A res t > 1548 A), as the blue half is affected by the mini- 
BAL troughs B and C. The Civ emission-line light curve 
shows some low-level variations throughout the SDSS- 
RM campaign, changing by about a factor of ~10% dur¬ 
ing this time period, but no large-scale changes similar 
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Rest-Frame Wavelength (A) 

1430 1440 1450 1460 1470 1480 1430 1440 1450 1460 1470 1480 



Fig. 9.— Four different pairs of continuum-normalized spectra between which the EW of Trough A varies at greater than 4a significance 
(A EW > 4 x ct^ew)- The continuum-normalized spectra are indicated in red and blue, and the horizontal dashed lines are at a normalized 
flux density of 1.0 and 0.9 to help guide the eye. The spectra have been smoothed for display purposes (only in this figure) using a boxcar 
window over 3 pixels to reduce noise. The vertical bars show the la uncertainties associated with each point in the smoothed, normalized 
spectra. 


Rest-Frame Wavelength (A) 

1430 1440 1450 1460 1470 1480 1430 1440 1450 1460 1470 1480 



Outflow Velocity (x 1000 km s l ) 


Fig. 10.— Flux ratios between the four different pairs of 
continuum-normalized spectra shown in Figure [9] The spectra 
have been binned by 3 pixels to reduce noise. Red vertical bars 
show the la uncertainty in the flux ratio in each binned pixel, 
and the horizontal bars show the size of the bins. Blue horizontal 
lines show the weighted mean of the flux ratio across the Trough A 
region. The x 2 °f the fit and the number of degrees of freedom 
(DOF) are shown in each panel. In all cases, a fit to this constant 
model is statistically acceptable. 



Fig. 11.— The decreasing part of the EW light curve for Troug h 
A (black filled circles) with a linear fit using the |Kelly| (2007) 
method that includes intrinsic scatter. The three pairs ot epochs 
between which we see significant variability are shown as colored 
squares (red, cyan, and yellow). The black solid line shows the 
linear fit, blue dotted lines show the measured intrinsic scatter, 
and red dashed lines show the expected scatter about the linear 
fit assuming no intrinsic scatter; i.e., if the expected scatter was 
only due to the uncertainties in EW. The expected scatter due to 
uncertainties in EW is smaller than that observed. 


to those shown by the BAL trough. However, the uncer¬ 
tainties are non-negligible due to the uncertainty of the 


spectrophotometric flux calibrations ( Shen et aT7||2015 ). 

Based off of its luminosity at rest-frame 1350 A (as 
measured from the mean spectrum) and the C IV R—L re- 
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Rest-frame Wavelength (A) 

1400 1450 1500 1550 



Outflow Velocity (km s *) 


Fig. 12.— Mean contin uum-normalized spectra in each group as 
defined in Section 13.2.31 of the text. Each color shows the mean 
of each group of spectra throughout the campaign. Horizontal 
dotted lines are normalized flux densities of 1.0 and 0.9 to guide 
the eye. The profile of Trough A remains the same throughout the 
campaign; only the strength/depth of the trough varies. 



Rest-frame days since MJD 56660.209 

Fig. 13.— The rest-frame EW of the Civ BAL (Trough A; top 
panel) and the integrated flux of the CIV emisson line (lower panel) 
as a function of rest-frame days since the beginning of the SDSS- 
RM campaign. The emission-line flux was calculated using only 
the red half of the Civ emission line, which is free of absorption 
features. 


lation of Kaspi et al. (2007), the expected mean timescale 
for C IV to respond to ionizing-flux changes is on the or¬ 
der of 130 rest-frame days, though gas nearer to the line 
of sight will respond more quickly. Because of the un¬ 
certainties in integrated emission-line fluxes and because 
the baseline of the current data from the SDSS-RM cam¬ 
paign is less than ~130 days, we are not able to explore 
fully the C iv emission-line variability — it might simply 
not have responded to the observed ionizing continuum 
yet. Additional observations of this target, with the ex¬ 
tended time baseline covered by the ongoing monitoring 
periods in 2015 and 2016, could be used to explore this 
further. 

We also consider the He n A1640 (n = 3 — 2) emission 
line as a potential probe of the ionizing continuum above 
54 eV. The Hen emission line is present but weak in our 
spectra; we m easured the Hell EW using the method of 
Baskin et al. (2013) for ease of comparison with their 
dataset. Our quasar’s average Hen EW of 3.5 ± 0.4 A is 
consistent with their BAL quasar average of 3.9 A, but 
lower than their non-BAL quasar average of 5.1 A (their 
Table 2). From the pre-campaign BOSS spectrum to the 
last SDSS-RM spectrum, Hen shows a marginal weak¬ 
ening in its rest-frame EW of approximately 20 ± 10%. 
The weakness of the He n emission line and the resulting 
large uncertainties in measuring its EW make it impossi¬ 
ble for us to measure short-term variability, and thus the 
He II emission in this target is not instructive as a probe 
of the short-term variability of the ionizing continuum. 

4. DISCUSSION 

4.1. Observed Variability in Context 

Our observations of C IV BAL variability in this quasar 
are the fastest reported detect ions of significan t BAL 
variability to date. Previously, Capellupo et al. (20131 
reported variability on timescales down to 8-10 days m 
two quasars; our shortest observed timescale is shorter 
by more than a f actor of five. The targets observed by 
Capellupo et al. (2013) showed variability in only small 


portions of their respective troughs; we here see coor¬ 
dinated variations across the entire velocity range of the 
trough, with no discernible BAL-profile changes through¬ 
out the campaign. We see very little significant variabil¬ 
ity in the rest-frame UV continuum as measured by the 
i-band photometric light curve and within the line-free 
regions of the spectra. 

To determine whether this particular quasar and the 
ClV BAL are atypical when compared to the general 
BAL population, we compare their vari ous properties 
with those of the BAL quasar sample of Gibson et al.| 
(2009). Our target has an absolute magnitude and red- 
shift that are fairly typical of the BAL population; how¬ 
ever, the quasar itself is less reddened than typical BAL 
quasars. The BAL trough itself is somewhat weak, with 
a smaller EW than many ClV BAL troughs, and the 
higher-velocity portion of the trough is deeper than the 
low-velocity portion - this shape is more common in 
weaker BALs. The velocity width of the trough, at 
4340 km s -1 , is not unusual, but is somewhat larger 
than the majority of CIV BAL troughs. The maximum 
velocity of the trough falls within the normal range of 
the broader BAL population, but the minimum velocity 
of the trough falls towards the high end of the distribu¬ 
tion (the trough is well-separated from the CIV emission¬ 
line). In summary, neither the quasar nor its BAL are 
strong outliers in any of the properties examined, though 
in some cases they do fall in less common regions of the 
distribution among BAL quasars. 

We also investigate whether the amplitude and 
timescale of variability seen during the SDSS-RM cam¬ 
paign is normal for a BAL quasar by comparing our ob¬ 
served variability with that of the large sample of|Filiz Ak 


et al. (2013). Figure 14 shows the change and fractional 


change in EW as a function of A f, reproduced from|Filiz] 


Ak et al. (2013) with our data overplotted in red. Our 
points lie within the envelope defined in that work, sug¬ 
gesting that the variability properties seen in our target 
may be common in quasars; our target is not an apparent 
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At (years) 

Fig. 14.— Change in EW (AEW; top panel) and fractional 
change in EW (AEW/<EW>; lower panel) as a fun ction of time 
between epochs (At), reproduced from Filiz Ak ct al. (2013) . Gray 
points show data from Filiz Ak ct al. (2013), Barlow 1993),|Lund-| 
|gren et al. | (|2007|) , and|tlibson et ahlJ2UUS|). Aed squares snow data 
tor each BLISS- AM epoch paired with every other epoch during the 
campaign. Values and uncertainties for the SDSS-RM data were 
calculated following the prescriptions defined in Section 3.3 ofIFilizI 
|Ak et al.| (201^. '-' 


outlier in its variability characteristics. The variability 
amplitudes on short timescales during our campaign are 
also consistent wit h those observed on sh ort timescales 
in a few targets by Capellupo et al. (20131. 

There are a number ol potential causes of BAL variabil¬ 
ity previously discussed in the literature, and the physi¬ 
cal constraints we can place on the nature of the outflow 
are dependent on the mechanism driving the variability. 
Below we briefly consider a few possible causes of the 
observed BAL variability and use our measurements to 
determine the relevant physical parameters of the out¬ 
flow. 


4.2. Possible Causes of Variability 
4.2.1. Ionization State Variability 

Recent studies by Filiz Ak et al. (2012, 2013) have 
demonstrated that, in general, distinct Civ BAL troughs 
in the same object separated by as much as 10,000-20,000 
km s -1 tend to strengthen or weaken together. Such cor¬ 
related changes over large velocity ranges suggest that 
much BAL variability arises from changes in the ionizing 
flux reaching the absorber rather than from transverse 
gas motions. There are several aspects of our observa¬ 
tions that suggest variations in the ionization state of 
the gas are responsible for the observed variability. First, 
global changes occur across the entire trough rather than 
in small segments — this can be understood if the ion¬ 
ization state of the gas changes but is difficult to explain 
in other scenarios (see below). Secondly, we see coor¬ 
dinated variability between the high-velocity ClV BAL 
Trough A and the mini-BAL Trough B (Figures [7] and 
[8|, which are separated by ~12,000 km s -1 . 

A notable aspect of our observations is the lack of sig¬ 
nificant variability in the rest-frame UV continuum flux, 
as measured in the photometric light curves (Figure [4]); 
at most, the photometric flux in our observations varies 
by only 10%. However, the flux in the i band is not 


necessarily a reliable tracer of the ionizing continuum 
C IV is ionized by radiation in the extreme ultraviolet 
(EUV) part of the spectrum, which is likely m ore vari¬ 
able t han the observed emission. For example, |MarshaIIl 
et al. ([1997) monitored the Seyfert 1 galaxy NGG 5548 
using the extreme Ultraviolet Explorer , finding that the 
EUV continuum at 0.15-0.2 keV was significantly more 
variable than that of the continuum at 1350 A; at times 
the EUV was observed to vary by factors of ~2 while the 
lower-energy UV as observed with the Hubble Space Tele¬ 
scope varied by only ~10% during the same period. The 
EUV was found to vary in phase with the UV/optical 
continuum to within the time sensitivity limit imposed 
by the cadence of the observations. 

While there are few other e xamples of AGN varia bility 
monitoring in the EUV (e.g., Marshall et al. 1996), the 
same trend has been seen in X-rays; simultaneous ob¬ 
servations in the X-ray and optical have shown that the 
short-term variations in X-ray flux are often not detected 
in the optical, t hough correlations ar e seen on longer 
timescales (e.g., 
were done with 


Peterson et al. 2000). These studies 
ocal Seyfert galaxies rather than dis¬ 
tant quasars; however, the X-ray variability properties 
are broadly similar among the two populations - both 
quasars and ne arby Seyferts show strong, rapid X-ray 
varia bility (e.g., Gibson & Brandt||2012 Shemmer et al.| 
2014|. Such strong variations of the ionizing continuum 
could be sufficient to drive the observed rapid ionization 
response of the BAL trough. Furthermore, if the ion¬ 
ization conditions of the absorber are such that the C IV 
abundance has a steep depen denc e on the i onization pa - 
rameter (e.g., see Figure 1 of Kallman & McCray|1982), 
even more modest changes in the ionizing continuum can 
drive strong changes in the C iv column density along the 
line of sight. 

Clump iness of the X-ray “shielding gas” (from the 


model of Murray et al. 


1995) is another possible expla¬ 


nation for the observed lack of photometric variability. 
In this model, the shielding gas is located between the 
continuum source and the outflow, and prohibits high- 
energy radiation (soft X-rays and potentially some of 
the ionizing continuum) from reaching the outflow, but 
does not necessarily affect the lower-energy observed UV 
continuum. If the shielding gas is corotating with the 
disk, dumpiness or movement of the shielding gas could 
cause the amount of ionizing flux reaching the outflow to 
vary on shorter timescales than the ionizing continuum 
changes themselves. Models do predict variabil ity on 
short timescales related to th e shielding gas (e.g., |Proga| 
et al. 2000 Sim et al. 2010) and a fe w studies report 
possible detections of this as well (e.g., [Gallagher et al.| 
2004 Saez et al.||2012). If our BAL variability is driven 
predominantly by changes in the shielding gas, then the 
rapidity of the trough variability constrains models of the 
size and kinematics of the shielding gas rather than the 
outflow gas. 

4.2.2. Cloud Crossing 

Another potential source of variability is the transverse 
movement of the outflow material across our line of sight 
to the quasar source (“cloud crossing”). Several previ¬ 
ous studies of BAL variability are consistent with this 
scenario, as variability is often o bserved in only narrow 
portions of the BAL troughs (e.g., Hall et al.|2011 Vivek 
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et al.|2012 Capellupo et al.|2013 1. However, for the C iv 
BAL trough in our target, which is ^4300 kms -1 wide, 
to undergo coherent strengthening or weakening across 
its entire velocity span as observed would require the ab¬ 
sorbing gas at all velocities to be physically connected. 
This explanation is implausible, as the trough width is 
too large to be sustained as a turbulent velocity disper¬ 
sion wi thout dispersing the a bsorbing structure (e.g., see 
§5.2 of Rogerson et al. |2011|) or heating it to a temper¬ 
ature at which no GIV would be seen (the thermal ve¬ 
locity in photoionized gas exhibiting C IV absorption will 
be < 10 kms -1 ). An absorbing structure with a gradi¬ 
ent in line-of-sight velocity across the continuum source 
can explain the trough width, but not the coherent vari¬ 
ability as a function of velocity, because gas at different 
velocities will cover different regions of the continuum 
source and will not be moving into or out of the same 
line of sight at the same time, as required to explain our 
observations. We can think of no plausible absorbing 
structure whose bulk motion, either into or out of our 
line of sight to the entire continuum emission region (or 
to hot spots within that region), can reproduce our ob¬ 
servations. Furthermore, this scenario does not explain 
the coordinated variability seen between Trough A and 
Trough B. 


4.2.3. Inhomogeneous Emission 

Partial coverage of a variable inhomogeneous accretion 
disk (e.g., Dexter & Ago! 2011) by an absorber can cause 
EW variability in the resulting trough. If the part of 
the disk not covered by an absorber brightens, both the 
observed continuum and the residual flux in the trough 
will increase, and the normalized trough will appear to 
weaken; conversely, if that part of the disk fades, the EW 
of the trough will increase due to the fading continuum. 
If the part of the disk covered by an absorber bright¬ 
ens, the observed continuum flux will increase more than 
the residual flux in the trough will, and the normalized 
trough will appear to strengthen; conversely, if that part 
of the disk fades, the normalized trough will appear to 
weaken. However, in any of the above cases, for the mean 
fractional depth of the trough to decrease smoothly by 
12% of the continuum level from rest-frame day 18 to 53, 
as observed, would require the observed continuum flux 
to change smoothly by a similar percentage in that time 
frame. The photometric light curves suggest no such 
trend. Therefore, we do not consider an inhomogeneous 
emitter a viable explanation for the observed variability. 


keeping C v fixed at all v, both models are acceptable fits 
to our data. However, the former scenario requires co¬ 
ordinated changes in C v across the entire width of the 
profile, for w hich we found no plausible mechanism in 
Section [4.2.2| On the other hand, the latter scenario is 
naturally explained by ionization variability. Adopting 
that scenario, we can constrain the gas to be at most 
marginally optically thick, with r max = 0.55 (for con¬ 
stant C v = 1) to r max = 1.5 (for constant C v = 0.47), 
and with gas at most velocities and in most epochs hav¬ 
ing t v well below r max . 


4.4. Physical Constraints from Ionization State 
Variability 

Following the discussion in Section |4.2[ we conclude 
that a change in ionization state of the absorbing gas is 
the most likely cause of the observed BAL variability in 
our target. Regardless of whether the change in ionizing 
flux reaching the absorber is caused by fluctuations in the 
high-energy ionizing continuum itself, or due to changes 
in shielding gas, we can use our observations to place 
constraints on the outflow properties. For gas which is 
initially in photoionization equilibrium, the characteris¬ 
tic timescale for changes in the density nt of ionization 
stage i of some element in response to an ionizing flux 
change can be written as 


t* = [~fi (Ii - n e oti-i)Y 


(1) 


(equivalent to Eq. 10 of Arav et al. 2012), where —1 < 
fi < Too is the fractional change in the ionization 
rate per ion of stage i ( Ift > 0) = (1 + fi)Ii{t = 0)), 
at -1 is the recombination coefficient to ionization stage 
i — 1 of the ion, and a negative timescale represents a 
decrease in ni with time. This equation only considers 
photoionization and recombination; other processes are 
neglected. Assuming an absorber of constant thickness, 
observed changes in ionic column density (i.e., r) corre¬ 
spond to ionic density changes. Gas showing such varia¬ 
tions is not in a steady state by definition, but such gas 
can still be in equilibrium with a varying ionizing flux 
if its t* is considerabl y sh o rter than the flux variability 
timescale (e.g., §6 of Krolik fc Kriss||1995 ). For optically 
thin gas at distance r from a quasar with luminosity 
at frequency v, the ionization rate per ion of stage i is 
given by 


4.3. Physical Constraints from the Lack of 
Trough-Profile Variability 

We found in Section [3. 2. 3| that Trough A is consistent 
with having an unchanged line profile throughout the 
campaign; only its strength varies. The normalized flux 
density at velocity v in a trough can be characterized 
as F v = (1 — C v ) T C v e ~ Tq ’, where C v is the frac¬ 
tion of the emission source covered by the absorber and 
t v is the optical depth in the absorbing ion, both along 
our line of sight. In cloud crossing variability, changes 
in C v predominate; in ionization variability, changes in 
t v predominate. The S/N of our spectra are not high 
enough to distinguish statistically between changing C v 
while keeping r„ fixed at all v and changing t v while 


where cr„ is the ionization cross-section for photons of 
energy hv. 

If the absorbing gas is sufficiently far from the 
quasar that /, -C n e at- i, then Equation 1 reduces to 
Gee = 1/ fiiT-eO-i-i (which is just the recombination 

time of the ion in the = — 1 case where the ioniz- 
ing flux drops to zero; e.g., Capellupo et al. 2013), and 
the observed absorption variability timescale constrains 
the density of the absorber. However, if the absorbing 
gas is close enough to the quasar that n e oti-\, then 
the relevant timescale is G = — 1 /fih and the absorption 
variations of the ion reflect the ionizing flux variations of 
the quasar, with no density constraint derivable just from 
absorption variations. An observed timescale for vari¬ 
ations in optically thin absorption therefore constrains 
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the absorbing gas to either have a density n e > n m i n 
and r > r eqU ai, where 7' e quai is the distance at which 
Ii = n m i n ai-i, or to be located at r < r equa i with no 
constraint on the density. 

To determine the constraints in our target, we assume a 
temperature of log T = 4.3, the tempe rature of t he broad 
line region gas in the Civ region (e.g., Krolik 1999). The 
recombination coefficient corresponding to this temper- 
ature is acm = 2.45 x 10~ n cm 3 s -1 (from the 

CHIANTI online database; |Dere et aL|l997j |Landi et al.| 
20131. For the simple case of the ionizing nux dropping 
to zero, fi = — 1 and the timescale t* above can be ap¬ 
proximated as the recombination time f rec ~ 1 /n e acin- 
Using our observed timescale of 1.20 days as an upper 
limit to t r ec for Trough A, we calculate a lower limit on 
the density of the gas to be n e xL 3.9 x 10 5 cm -3 . We 
cannot search Trough B for short-term variations as we 
can in Trough A, because Trough B sits on top of a slowly 
varying, imperfectly modeled emission line. However, we 
can set a lower limit on n e for Trough B using the ob¬ 
served EW decrease of about 20% in 20 days. This indi¬ 
cates a characteristic timescale of ^100 days and a min¬ 
imum n e about 100 times lower than that in the Trough 
A absorber. 

Our calculation for Trough A at log T = 4.3 yields sim- 
ilar densities to those found by Capellupo et al. (20131 
for their targets, which showed variability on timescales 
down to 8-10 days, though they used a different recombi¬ 
nation coefficient (presumably calculated for a different 
temperature). Using the same recombination coefficient 
they used with our observed timescale would yield n e for 
our target higher than their measurement by more than 
a factor of six [n e 2.7 x 10 6 cm -3 ). 

Using our lower limit density of n e (J, 3.9 x 10 5 cm -3 , 
we calculate the minimum distance from the quasar at 
which that lower limit is valid. From its observed flux 
density at rest-frame 3000 A, our quasar has a bolometric 
luminosity AboI = 1.7 x 10 47 erg s -1 using a bolometric 
correction factor of 5 (Richards et al. 20061 . We adopt 
the U V-soft spectral energy distribution of Dunn et al. 
(2010) to calculate and integrate Equation[2| setting 
U = n e a cm to solve for the radius ?' eq uai at whicn the ion¬ 
ization rate is equal to the recombination rate. We find 
equal = 3.8x 10 2 ° cm (120 pc). Thus, the absorber in this 
quasar either has n e (i, 3.9 x 10 5 cm -3 and r > 120 pc, 
or it has r < 120 pc, in which case the recombination 
timescale cannot be used to constrain the density. This 
radius is much larger than the launching radius proposed 


outflows at kpc-scale distances, which are on the order of 


(e-fG 

Murray et al. 1995 

Proga et al. 2000’ 

~-— ~ 

Higginbottom 

et al. 

2013) — some also theorize that BA1 

jS are launched 


emitting region (e.g., Baskin et al. 2014). A physical 


mechanism for launching such high-velocity outflows at 
distances of parsecs or kiloparsecs that is consistent with 
observations is difficult to produce. However, r > 120 pc 
falls within the wide spread of outflow radii reported in 
the literature — many studies report outflow radii on 
the o rde r of parse cs or kiloparsecs (see, e.g., Table 10 of 


Dunn et al. 2010 and references therein), thou gh some 
studies do report significantly sma ller radii (e.g., Leighly 
et al.||2009 Capellupo et al.||2013 ). Our measured lower 
limit on n e is higher than those typically found for the 


10 3 — 10 4 cm 3 (e.g ., Dunn et al.||2010 Hall et aJ7||201 1 


Borguet et al. 2013), though we make this comparison 
with cautloiTaue to the differences in analysis methods 
and possible biases within the samples of outflows. Our 
density lower limit is closer to other measurements of out¬ 
flows that are found at distances on the order of r equa i 
(~100 pc), which have n e _ ~ 10 4 — 10 5 cm -3 (e.g., Arav| 
et al.|2013) Borguet et al.||2013 ). 


The discrepancy between observations and model pre¬ 
dictions can likely be reconciled through selection effects 
that lead to at least some BAL outflows (particularly 
those with distance measurements to date) being ob¬ 
served at ra dii much larger than the launching radii (e.g., 
Lucy et al. ~ 20141). The radius at which BAL outflows are 


observed is not necessarily the radius from which they 
are launched; some studies suggest that outflows can 
maintain high velocities at large distances allowing us 
to observe t hem much farther out than the launching lo¬ 
cation (e.g., Faucher-Giguere et al. 2012). In fact, we can 
use our observations to estimate tne upper limit for the 
launching radius Haunch of the outflow; however, to do so 
requires the adoption of an outflow model. For example, 
within the conte xt of the radiat ion pressure confinement 
(RPC) model of Baskin et al. (2014), the terminal ve¬ 


locity of an outflow is n te rm ~ 7'Kep x/rA/LEdd , where r 
is the force multiplier and nxep is the circular velocity 
of gas at the launch ra dius (e.g., Laor fe B randt 2002). 
Adopting r ~ 10 from Baskin et al. "(20141) and assum- 

ii 


10 from 

ing that our BAL is observed at terminal velocity, we 


find ri aunc h < 0.3 pc ~ 400 light-days. The mean ra¬ 
dius of the C iv-emitting broad line region in this quasar 
is est imated to be Rblr ~ 130 liglrt-days (see Section 


3.2.4), and with Mbh = 1-1 x 10 9 Mg, this target has 
a Schwarzcliild radius -Rsch ~ 0.125 light-da ys. Thus for 
this quasar, assuming the RPC model of |Baskin et al.| 
(|2014 1 yields Haunch < 3I? B lr ~ 3200-Rsch- Note that 


typically 


Haunch is model dependent and thus different values will 
be obtained depending upon the model. 

We can also place constraints on n e using the pres¬ 
ence of other absorption features, or the lack thereof. 
For example: With a density of n e xf 3.9 x 10 5 cm -3 
at r = r equa i, the gas would have an ionization 
parameter of Uh = 0.07 for our assumed SED 
(Uh = Qh/^^uhC, with Qh = l-2x 10 57 hydrogen- 
ionizing photons s -1 and uh = 0.82n e ). That value of 
Uh is consistent with the presence of absorption from 
Civ but not CiiA 1335 in this object. A factor of ~50 
lower Uh would yield CiiA 1335 absorption about hal f 
as deep as that of Civ (Fig. 3 of Hamann et al. 1995), 
which is not observed in this source. Given that Uh 
can be lower by at most a factor of 50, if the gas is lo¬ 
cated at r > r equa i, it is thus also constrained to have 
7’equai < r < x/507-equai and also must have density 
3.9 x 10 s <n e < 10 7 (r equal /r) 2 cm -3 . 

We next measure the kinetic luminosity of the outflow 
relative to its Eddington luminosity to evaluate the po¬ 
tential of the outflow to provide substantial feedback to 
the quasar host galaxy. To estimate the Eddington lumi¬ 
nosity (Le<m) of this quasar, we first require an estimate 
of Mbh- We use Equation 9 from Rafiee & Hall (2011), 
which incorporates the dispersion er (in km s -i ) of the 
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Mg II line: 


14410 . 5^2 


Mbh/M© = 30.5[AL3ooo/10 

FWHM/2.35 


We assume a = 

FWHM = 3274 km s" 1 


(Paris et al. 2015). 


(3) 

with 

The 


observed flux density at rest-frame wavelength 3000 A is 
/3000 = 8xl0 -17 erg s -1 cm -2 A~ 4 , which corresponds 
to AL 3000 = 3.38 x 10 46 erg s _1 . Using the above 
equation yields Mrt t = 1.1 x 1 0° Mm with an intrinsic 
scatter of 0.4 x 10 9 ( Rafiee fc Hall|[20lT I. With our cal¬ 
culation of Mbh above, we measure LEdd = 1.7 x 10 47 
erg s _1 , and thus in this case the Eddington ratio, 
/Edd = iboi/^Edd, is about 1. We note that Mbh 
measurements made from single spectra such as this are 
subject to a number of uncertainties; we consider our 
measured Lecm to be an approximation. 

The kinetic luminosity of a BAL outflow consisting of 
a thin shell at distance r from the quasar is 

Lk = 2TTfj,m p QrNHV 3 (4) 

(e.g., Eq. 1 of Arav et al. 2013) where fl is the global 
BAL covering fraction (unity for a spherical outflow), 
fj,m p = 1.4m p is the mean mass per proton in the ab¬ 
sorber, and Njj is the absor ber’s hydrogen column. We 
assume C v = 1 (see Section 14.31) and H = 0.2 (derived 
fro m the overall fractio n of quasars hosti ng Civ BALs; 
e.g., Gibson et al.|2009 or Dunn et al. 2010), and the aver¬ 


age (J IV t v of Trough A is 0.26. The column density per 
km s _1 for an unsaturated and resolved doublet a bsorp¬ 
tion feature is N v = (3.7679 x 10 14 cm _2 )T„/A/,j (Arav 
et al. 1999 Hall et al. 2003); we assume an unresolved, 
optically thin trough and use an average A = 1549.055 A 
and a summed oscillator strength f l3 = 0.286. For 
Trough A, with a total width of 4340 km s _1 , we have a 
total Acrv = 9.6x 10 14 cm~ 2 . Conservatively assuming 



Lk/LEdd = 1-1 x 10 4 , or 0.01%. AGN/host-galaxy 
feedback models require Lk/LEdd A 0.5% to 5% in order 
for the outflowing material to provide substantial feed¬ 


back to the host galaxy (e.g.,|Di Matte o et al. 2005 Hop¬ 
kins et al. 2005 Hopkins & Elvis 2010), which is about an 


order of magnitude larger than our measurement. Thus, 
we find no evidence that the outflow will contribute sig¬ 
nificant feedback to the host galaxy. However, it could 
do so if it is accompanied by a sufficient column density 
of very high ionization gas, to which we are not sensitive. 


5. SUMMARY AND FUTURE WORK 

We have discovered and studied a rapidly variable C IV 
BAL trough in a quasar observed as a part of the SDSS- 
RM campaign using 31 epochs of BOSS spectroscopic 
observations and accompanying optical photometric light 
curves obtained over six months in 2014. Our main find¬ 
ings are the following: 

1. We measure significant variability of the EW of 
the high-velocity Civ BAL trough (Trough A) on 
timescales as short as 1.20 days in the quasar rest 


frame, the shortest timescale yet reported. The 
EW varied on the order of ~10% on the smallest 
timescales, while over the entire campaign, the EW 
varied by more than factor of two. See Section [3~2j 

2. The rest-frame UV photometric flux of the quasar 
shows little variability (~10%) during the length of 
the campaign. See Section \2A\ 

3. The trough itself shows coordinated variability 
across its entire velocity range throughout the cam¬ 
paign, with no significant change in shape across its 
4340 km s _1 span, suggesting that the variability 
is not due to bulk motion of the gas. See Sections 
1372731 and l4T2~2l 

4. We measure coordinated variability between a 
lower-velocity C IV mini-BAL (Trough B) and the 
highest-velocity ClV BAL (Trough A); this be¬ 
havior, combined with the coordinated variability 
across the entire Trough A structure, suggests that 
the variability is driven by a chang e i n ioniz ation 
state of the gas. See Sections |3.2.1| and |4.2.1| 

5. For a change in ionization state, a variability 
timescale of 1.20 days suggests a lower limit to 
the density of the gas of n e ^ 3.9 x 10 s cm -3 
so long as the gas is at a radius of greater than 
~ 120 pc. We also measure the kinetic luminos¬ 
ity Lk = 2.2 x 10 -4 LEdd, which suggests that 
the outflow itself is not a major contributor to host 
galaxy feedback; however, since this is an approxi¬ 
mation only , we do not rule out the possibility. See 
Section 14.41 


The key to this successful investigation of such short- 
timescale variability was the high cadence of the SDSS- 
RM spectroscopic data; this allowed a robust exploration 
of changes on shorter timescales than have previously 
been investigated. There are several other BAL quasars 
that were observed as a part of this program that have 
not yet been thoroughly investigated; a search for sim¬ 
ilar short-term BAL variability in these targets will be 
instructive and will help determine the frequency of this 
behavior. Additional monitoring of these BAL quasars 
will be obtained in the following years as a part of the 
ongoing SDSS-RM program. While the upcoming obser¬ 
vations will not provide spectra at the same high cadence 
that was obtained in 2014, the additional data will allow 
a search for continued BAL variability in these targets 
and extend the baseline for the longer-term trends ob¬ 
served (for example, the troughs in this object are present 
at approximately the same strength in the first few spec¬ 
tra of 2015). In addition, the extended baseline of these 
observations will allow us to explore the relationship be- 
twee n BAL trough and emission-line variability (see Sec¬ 


tion 

due 


3.2.4), which has not previously been investigated 
o the lack of suitable data for such a study. 

Figure [TT| shows that the variability properties of this 
target are consistent with those found in other quasars, 
indicating that short-term variability such as we see may 
be common. High-cadence, high-S/N spectroscopic ob¬ 
serving programs geared toward investigating variability 
in additional BAL quasars would expand the sample, al¬ 
low us to constrain further the dominant driving mech¬ 
anisms of BAL variability in quasars, and refine models 
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describing the environment and physics regulating BALs 
as well as their possible contributions to feedback to their 
host galaxies. 
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